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Abstract 



Nuclear-size corrections of order (Za) 5 and (Zo) 6 to the S-state levels 
of hydrogenic atoms are considered. These nuclear-elastic contributions are 
somewhat smaller than the polarizability (nuclear-inelastic) corrections for deu- 
terium, but are of comparable or larger size for the hydrogen case. For deu- 
terium the (attractive) nonrelativistic [Zaf 1 correction to the 2S-1S transition 
is 0.49 kHz, while the (repulsive) relativistic (Za) e contribution is —3.40 kHz. 
For the proton the corresponding corrections are 0.03 kHz and —0.61 kHz, re- 
spectively. The (Za) 5 contribution largely cancels the Coulomb-retardation 
part of the nuclear-polarization correction. 
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Introduction 



Recent experiments || |3| have pushed the precision of the spectroscopy of the 
hydrogen isotopes to the point where previously unconsidered higher-order terms (in 
a, the fine-structure constant) are now required. In the isotope shift of the 2S- 
1S transition between deuterium and normal hydrogen, for example, experimental 
uncertainties]^] of less than 0.2 kHz (one part in 3 • 10 4 of the total nuclear-finite- 
size correction) have been reported. The precision is such that these measurements 
afford us an unparalleled opportunity to extract a precise value of the d-p mean- 
square-radius difference. At the reported level of experimental precision, this quantity 
is sensitive to "exotic" contributions to the charge density^, such as relativistic 
corrections and meson-exchange currents. It therefore behooves us to calculate all 
higher-order contributions to the (atomic) frequency shift of this size. 

We report here a calculation of nuclear finite-size (i.e., nuclear-elastic) corrections 
of orders (Za) 5 and (Za) 6 , which supplement the usual leading-order (Za) A term. 
Recently, the polarization (nuclear-inelastic) corrections of order (Za) 5 were calcu- 
lated for deuterium by expanding that quantity in a series in ER ~ where E is 
the average (virtual) nuclear-excitation energy and R is a typical nuclear size. The 
first three orders in this series are expected to be accurate at the level of roughly 
0.01 kHz, although individual nuclear observables in that expression cannot be deter- 
mined to that accuracy]?]]. Very recently, the leading-logarithm contributions to the 
proton-polarizability correction were calculated ||, as well. 

Fortunately, little analytic work is required to obtain the finite-size corrections, 
since they were calculated many years ago in the context of muonic atoms || fTOfl . The 
first three orders of corrections for the nth S-state can be written in the form: 

AE n = ^Za |0 n (O)| 2 (V) - ^V> (2) + (Za) 2 F REL + (Zap) 2 F NR + ■ ■ 

where Z is the nuclear charge, (r m ) is the mth moment of the nuclear charge distri- 
bution (normalized to unit charge), p is the reduced mass, (fi n (0) is the electron wave 
function at the origin, and the Zemach moment|TT, |J (r 3 )(2) is defined by 



(r m ) (2) = J d 3 rr m p (2) (r), (2) 
where the convoluted (Zemach) charge density is given by 

P( 2 )(r)= f d 3 zp(\z-r\)p(z) = p®p. (3) 
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The nonrelativistic correction Fnr (of relative order (Za) 2 m 2 R 2 , where m e is the 
electron mass) is negligible and will not be considered further, while the corresponding 
relativistic correction is defined by 



REL 



-(r 2 ) <m(/3r)> + 



■0(n) + 7 - 



(5n + 9)(n- 1) 
4n 2 



+ 7 - 2 + /] 



REL 



(4) 



where ip(n) is the digamma function, 7 is Euler's constant, (3 = 2Zafi/n, and 
(r 3 )(l/r) 



'REL 



+ / rf 3 Sp(s) / 0? 3 tp(t) 0(s-t) 



y-3 „3 „2 _ f2 

(t 2 + s 2 ) ln(t/s) -— + - + - 



+ 6 / d 3 up(u) / dHp(t) / rf 3 sp(s) 9(u - t) Q(t - s) 



+ V7+-H 



s^t 2 2s 2 t 



9 v t w y 36m 2 9m 9 



(5) 



The n-dependent terms in Eq. (4) were calculated independently by Karshenboim[|12 



We first treat -Fr EL before discussing the smaller, nonrelativistic [Zotf contribu- 
tion. We calculate (r 2 ), (ln(2a/ir)), and Jrel using the same techniques adopted 
earlier to treat the nuclear-polarization observables||. The proton charge distribu- 
tion is taken to be an exponential (dipole form factor) with a radiusfO of 0.862 fm. 
All integrals are known analytically for this caseQ. The deuteron charge distribution 



is obtained by first solving for the deuteron wave function using a variety|14|, [15], [16], 
I7| , [IB] , ^D) of first-generation (i.e., older) and second- generation |21], ^3| (i.e., 
recent) potential models. The latter fit the nucleon-nucleon scattering data very well; 
the best of them can be considered as alternative phase-shift analyses. This "bare- 
deuteron" density is folded with the sum of proton and neutron densities, and a spline 
fit is performed on the result. The neutron density is taken to be that of a dipole form 
factor multiplied by Ag 2 , with A adjusted to fit the observed ||24|| neutron mean-square 
charge radius: —0.114 fm 2 . Finally, double and triple integrals are performed using 
the spline-fitted folded density (viz., the p's in Eq. (5)). 

The results are given in Table 1. We summarize the second-generation results for 
the deuteron as 



-^REL 

(ln(2a/ir)) 
z/ (6) 



-3.094(4) fm 2 
-9.773(3) 
-3.40(1) kHz. 



(6) 
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Table 1: Deuteron finite-size corrections of order {Za) e for various potential models. 
The mean-square radius of each potential model is (r 2 ), the logarithmic radius is 
defined by (ln(2a/ir)), the relativistic density correlation is labelled Jrel, and the 
corresponding deuteron 2S-1S finite-size frequency shift is z/ F 6 g. The corresponding 
proton case is considered last using an exponential charge distribution. Both shifts 
are repulsive. 



Potential Model 


(r 2 ) (fm 2 ) 


(ln(2a/xr)) 


Jrel (fm 2 ) 


4 6 i (kHz) 


Second-Generation Potentials 


Argonne Vis 
Reid Soft Core (93) 
Nijmegen (loc-rel) 
Nijmegen (loc-nr) 
Nijmegen (nl-rel) 
Nijmegen (nl-nr) 
Nijmegen (full-rel) 


4.507 
4.505 
4.507 
4.498 
4.496 
4.494 
4.483 


-9.770 
-9.771 
-9.771 
-9.773 
-9.773 
-9.776 
-9.771 


-3.091 
-3.095 
-3.091 
-3.090 
-3.096 
-3.094 
-3.098 


-3.41 
-3.41 
-3.41 
-3.41 
-3.40 
-3.40 
-3.39 


First-Generation 


Potentials 






Bonn (CS) 

Argonne V14 

Nijmegen (78) 

Super Soft Core (C) 

de Tourreil-Rouben-Sprung 

Paris 

Reid Soft Core (68) 


4.505 
4.556 
4.579 
4.595 
4.530 
4.516 
4.459 


-9.772 
-9.762 
-9.757 
-9.755 
-9.766 
-9.768 
-9.774 


-3.099 
-3.108 
-3.111 
-3.116 
-3.098 
-3.094 
-3.056 


-3.41 
-3.45 
-3.46 
-3.48 
-3.43 
-3.42 
-3.38 


Proton 


proton 


0.743 


-10.652 


-0.473 


-0.61 



We caution that these "uncertainties" are merely spreads in the potential-model re- 
sults. Changing various aspects of the models that we have used (including nucleon 
radii) could produce larger changes than these uncertainties. 

As an example of this caveat we note that the small Darwin-Foldy relativistic 
correction to the proton and deuteron charge densities have not been included in 
our model. This is easily done by adding 0.0332 fm 2 to the (r 2 ) of the proton (and 
hence to the deuteron). Consequently, the effect cancels H] in the isotopic difference. 
This addition modifies the proton charge radius to 0.881 fm, which increases i/ FS by 
0.03 kHz in both cases, largely by increasing (r 2 ) rather than by changing Jrel or 
(ln(2a//r)). 
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Finally, we consider the (Za) 5 (second) term in Eq. (1). This is very similar in 
structure to the Coulomb-retardation part of the nuclear-polarization correction |J, 
which can be written in the form 



7T 2 

AE,„ = a m P 



3 



J d?x J rf 3 2/ |x-y| 3 (0|pt(x)p(y)|0)-Z 2 (r 3 ) 



(2) 



(7) 



where the functions p(x) and p(y) in the correlation function (0|pt(x)p(y)|0) are 
nuclear charge operators. Ignoring the difference between m e and p (which is of recoil 
order), the last term in Eq. (7) exactly cancels the second term in Eq. (1), leaving only 
the correlation function. That is, the separation of the nuclear Compton amplitude 
(which determines many of the nuclear-structure-dependent atomic corrections) into 
nuclear-elastic and nuclear-inelastic parts is somewhat artificial for this particular 
term. 



The correlation function in Eq. (7) vanishes for certain (special) cases of interest 
to us. For the deuteron case the correlation function generates two small terms from 
the finite sizes of the proton (2.44 fm 3 — > 0.032 kHz) and the neutron (—1.58 fm 3 — > 
—0.020 kHz), which largely cancel and leave a very small residue: 0.012 kHz. Both 
terms vanish in the point-nucleon limit because (Y^i^i) 2 = J2i^ 2 , where counts 
the charge of the ith nucleon. We note that the same cancellation takes place if 
we evaluate the proton correlation function in the naive, nonrelativistic, pointlike 
quark model (where mutatis mutandis counts the quark charges). In this case the 
individual cancelling quantities are only 0.03 kHz in size. 



Table 2: Contributions in kHz to the higher-order deuteron and proton finite-size 
frequency shifts for the 2S-1S transition, together with their order in (Za), differ- 
ences between the deuteron and proton, and grand totals. Negative contributions are 
repulsive. For comparison the Coulomb-retardation nuclear-polarization correction, 
f^j, is listed also, but is not included in totals. 



Order 


ret 
U pol 


(Zaf (Za) 6 


Total 


deuteron 
proton 


-0.48 
-0.03 


0.49 -3.40 
0.03 -0.61 


-2.91 
-0.58 


Isotope Shift 




0.46 -2.79 


-2.33 



These results are displayed in Table 2, where the finite-size Coulomb corrections 
are listed on the right, while the Coulomb-retardation contributions to the nuclear- 
polarization correction are listed (for comparison only) on the left (they are not 
included in the totals). The cancellations are evident. 
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In summary, we have computed the Coulomb nuclear-finite-size corrections of or- 
ders (Za) 5 and (Za) 6 to the energy levels of hydrogenic atoms. The relatively small 
(Za) 5 term largely cancels the previously-calculated Coulomb-retardation term to 
the nuclear-polarization corrections for the deuteron and proton. The {Za) 6 contri- 
bution is larger because of a very large logarithm and because its relativistic origin 
precludes the need for the very small nonrelativistic factors of (m e R ~ 255). The net 
higher-order Coulomb nuclear- finite-size corrections are —2.91 kHz for deuterium and 
—0.58 kHz for the proton. 
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